December 17, 2008 1:5 



WSPC - Proceedings Trim Size: 9in x 6in icap2008-Biraben 



00 

O 
O 

<N 

o 

Q 

\Q , Determination of the fine structure constant with atom 

interferometry and Bloch oscillations 

< — ] ] M. CADORET, E. de MIRANDES, P. CLADE, C. SCHWOB, F. NEZ, L. JULIEN 

and F. BIRABEN 

Laboratoire Kastler Brossel, ENS, CNRS, UPMC, 4 place Jussieu, 
75252 Paris CEDEX 05, France 
* E-mail: biraben@spectro.jussieu.fr 
' www.spectro.jussieu.fr 

Q ' S. GUELLATI-KHELIFA 

• Institut National de Metrologie, Conservatoire National des Arts et Metiers, 

^"V 61 rue Landy, 93210 La Plaine Saint Denis, France 

We use Bloch oscillations to transfer coherently many photon momenta to 
atoms. Then we can measure accurately the recoil velocity hk/m and deduce 
the fine structure constant a. The velocity variation due to Bloch oscillations 
is measured using atom interferometry. This method yields to a value of the 
0^ , fine structure constant a -1 = 137.035 999 45(62) with a relative uncertainty 

. of about 4.5 x 1CT 9 . 
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1. Introduction 

The fine structure constant a is the fundamental physical constant charac- 
terizing the strength of the electromagnetic interaction. It is a dimensionless 
quantity, i.e. independent of the system of units used. It is defined as: 

„2 

(1) 



47r eahc 



where eo lS the permittivity of vacuum, c is the speed of light, e is the 
electron charge and h = h/2ir is the reduced Planck constant. The fine 
structure constant is a key of the adjustment of the fundamental physical 
constants. 1,2 The different measurements of a are shown on Fig. 1. These 
values are obtained from experiments in different domains of physics, as the 
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quantum Hall effect and Joscphson effect in solid state physics, or the mea- 
surement of the muonium hyperfine structure in atomic physics. The most 
precise determinations of a are deduced from the measurements of the elec- 
tron anomaly a e made in the eighties at the University of Washington 3 and, 
recently, at Harvard. 4-6 This last experiment and an impressive improve- 
ment of QED calculations 7 ' 8 have lead to a new determination of a with a 
relative uncertainty of 3.7 x 1CP 10 . Nevertheless this last determination of 
a relies on very difficult QED calculations. To test it, other determinations 
of a are required, such as the values deduced from the measurements of 
h/mcs 9 and ft./ m Rb {mcs and mRb are the mass of Cesium and Rubidium 
atoms) which are also indicated on the figure 1. In this paper, we present 
the measurements of h/m^b made in Paris in 2005 and 2008. 

The principle of our experiment is the measurement of the recoil velocity 
v r of a Rubidium atom absorbing or emitting a photon (v r = hk/m , where 
k is the wave vector of the photon absorbed by the atom of mass m). As the 
relative atomic masses A r are measured very precisely, the measurement of 
h/rrifib is a way to determine accurately a via the Rydberg constant .Roc,: 

2 2R oa A r (Rb) h 

a = T~r~\ v 2 ) 

c A r (e) TO Rb 

In this equation, the relative atomic mass of the electron A r {e) and the 
Rubidium A r (Rb) are known with the relative uncertainties of 4.4 x 10~ 10 
and 2.0 x 10~ 10 respectively. 10,11 As the fractional uncertainty of Roo is 
7 x 10~ 12 , 12,13 the factor limiting the accuracy of a is the ratio h/m^b- 

2. Principle of the experiment 

The principle of the experiment is to coherently transfer as many recoils as 
possible to the atoms (i.e. to accelerate them) and to measure the final ve- 
locity distribution. In our experiment, the atoms are efficiently accelerated 
by means of N Bloch oscillations (BO). The velocity selection and velocity 
measurement are done with Raman transitions. The experiment develops in 
three steps, i) Firstly, we select from a cold atomic cloud of 87 Rb a bunch of 
atoms with a very narrow velocity distribution. This selection is performed 
by a Doppler velocity sensitive counter-propagating Raman transition. In 
2005, we have used a 7r-pulsc to transfer the atoms from the F = 2 to the 
F — 1 hyperfine levels of 87 Rb. In 2008, we have modified the experimental 
scheme to take advantage of Ramsey spectroscopy: we use a pair of tt/2 
pulses which produces a fringe pattern in the velocity space, ii) Secondly, 
we transfer to these selected atoms as many recoils as possible by means of 
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Fig. 1. Determinations of the fine structure constant in different domains of physics. 
The most precise measurements are shown in the lower part of the figure. They are 
deduced from the anomaly of electron and from the ratio h/mcs and h/m-^. We have 
taken into account the most recent result of the group of Gabrielse. 6 The two values 
deduced of h/m^ are presented in this paper. 



Bloch oscillations as explained later, iii) Finally, we measure the final veloc- 
ity of the atoms by a second Raman transition which transfers the atoms 
from the F = 1 to the F = 2 hyperfine level. In short, we have used two 
different pulse sequences, the it — BO — tt and ir/2 — ir/2 — BO — tt/2 — tt/2 
configurations. 

The accuracy of our measurement of the recoil velocity relies in the 
number of recoils (2JV) that we are able to transfer to the atoms. Indeed, if 
we measure the final velocity with an accuracy of a v , the accuracy on the 
recoil velocity measurement a Vr is: 

Bloch oscillations have been first observed in atomic physics by the 
groups of Salomon in Paris and Raizen in Austin. 14-16 In a simple way, 
Bloch oscillations can be seen as Raman transitions where the atom begins 
and ends in the same energy level, so that its internal state (F = 1 for 
87 Rb) is unchanged while its velocity has increased by 2v r per Bloch oscil- 
lation. This is illustrated on Fig. 2 which shows the atomic kinetic energy 
versus the atomic momentum. The velocity distribution obtained after the 
7r/2 — 7r/2 selection is also represented. Bloch oscillations are produced in 
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Fig. 2. Acceleration of cold atoms with a frequency chirped standing wave. The vari- 
ation of energy versus momentum in the laboratory frame is given by a parabola. The 
energy of the atoms increases by the quantity 4(2j + l)E r at each cycle. The Ramsey 
fringe patterns represents the momenta distribution of the atoms in the F = 1 hypcrfinc 
level. 

a one dimensional optical lattice which is accelerated by linearly sweeping 
the relative frequency of two counter propagating laser beams (frequencies 
v\ and v-i). The frequency difference Av is increased so that, because of the 
Doppler effect, the beams are periodically resonant with the same atoms 
{Av = 4{2j + l)E r /h, j = 0,1,2,3.. where E r Jh is the recoil energy in 
frequency units and j the number of transitions). This leads to a succession 
of rapid adiabatic passages between momentum states differing by 2hk. In 
the solid-state physics approach, this phenomenon is known as Bloch oscil- 
lations in the fundamental energy band of a periodic optical potential. The 
atoms are subject to a constant inertial force obtained by the introduction 
of the tunable frequency difference Av between the two waves that create 
the optical potential. 14 

We now describe the acceleration process following the Bloch formalism. 
If, after the selection, the atom has a well defined momentum hqo with 
| go I < k, the atomic wave function is modified when the optical potential 
is increased adiabatically (without acceleration) and becomes in the first 
energy band: 



/ 

with I G Z. Here |go) designs the ket associated to a plane wave of mo- 
mentum go and the amplitudes (fio correspond to the Wannier function 17 




(4) 
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in momentum space of the first band. When the potential depth is close 
to zero, the limit of the Wannier function <f>Q is 1 over the first Brillouin 
zone and zero outside. On the contrary, if the potential depth is large, the 
Wannier function selects several components in the velocity space. When 
the optical lattice is accelerated adiabatically, the Wannier function is con- 
tinuously shifted in the momentum space following the relation: 

|*(t)> = J2 Mlo + 2lk - mv(t)/h)\q + 2lk) (5) 
l 

where v(t) is the velocity of the optical lattice. The enveloping Wannier 
function <j>$ is shifted by mv(t) in momentum space. After the accelera- 
tion, the potential depth is decreased adiabatically and, in equation 5, the 
Wannier function selects only one component of the velocity distribution. 
At the end, the wave function is = \qo + 2Nk). If Av is the velocity 
variation due to the acceleration, the number of Bloch oscillations N is 
such as \hqo + mAv — 2Nhk\ < hk. Consequently, if the initial atomic ve- 
locity distribution fits the first Brillouin zone, it is exactly shifted by 2Nv r 
without deformation, as it is shown, on Fig. 2, for the velocity distribution 
produced by a pair of it/ 2 pulses. 

3. Results in the tv — BO — 7r configuration 

Our experimental setup has been previously described in detail. 18 ' 19 Briefly, 
we use a magneto-optical trap (MOT) and an optical molasses to cool 
the atoms to about 3 [iK. The determination of the velocity distribu- 
tion is performed using a ir — ir pulses sequence of two vertical counter- 
propagating laser beams (Raman beams): 20 the first pulse with a fixed 
frequency S se i, transfers atoms from 5Si/ 2 , \F = 2, rap = 0) state to 5Si/ 2 , 
\F = 1, m? = 0) state, into a narrow velocity class (width of about v r /15). 
Then a laser beam resonant with the 55*1/2 (F = 2) to 5Pi/ 2 (F = 3) 
cycling transition pushes away the atoms remaining in the ground state 
F = 2. Atoms in the state F = 1 make N Bloch oscillations in a vertical 
accelerated optical lattice. We then perform the final velocity measurement 
using the second Raman 7r-pulse, whose frequency is <5 meas . The popula- 
tions of the F = 1 and F = 2 levels are measured separately by using a 
time of flight technique. To plot the final velocity distribution we repeat 
this procedure by scanning tlic Rciniaii bccim frequency S meas of the second 
pulse. 

To avoid spontaneous emission and to reduce other stray effects (light 
shifts and refraction index), the Raman lasers and the optical lattice are 
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Fig. 3. Velocity spectra obtained in the n — BO — tt configuration. Here TVi and JV2 are 
respectively the number of atoms in F = f and F = 2 after the acceleration process. 
These two spectra are obtained by performing the Bloch acceleration upwards or down- 
wards. The frequency difference between these spectra corresponds 1780 recoil velocities. 



blue detuned by ~ 1 THz and ~ 40 GHz respectively from the one photon 
transition. The delay between the two 7r-pulses is 12 ms and their dura- 
tion 3.4 ms. The optical potential depth is 70 E r . For an acceleration of 
about 2000 ms -2 we transfer about 900 recoil momenta in 3 ms with an 
efficiency greater than 99.97% per recoil. To avoid atoms from reaching 
the upper windows of the vacuum chamber, we use a double acceleration 
scheme: instead of selecting atoms at rest, we first accelerate them using 
Bloch oscillations and then we make the three steps sequence: selection- 
acceleration- measurement. This way, the atomic velocity at the measure- 
ment step is close to zero. In order to eliminate the effect of gravity, we 
make a differential measurement by accelerating the atoms in opposite di- 
rections (up and down trajectories) keeping the same delay between the 
selection and measurement 7r-pulses. The ratio h/m can then be deduced 
from the formula: 



m 



. _/v do 



down. 



(6) 



2(N u p + N down )k B (ki + k 2 ) 

where (S meas — 5 se i) up / down corresponds respectively to the center of the 
final velocity distribution for the up and the down trajectories, j\[ u p/ down 
are the number of Bloch oscillations in both opposite directions, fcs is the 
Bloch wave vector and k\ and k 2 are the wave vectors of the Raman beams. 
Moreover, the contribution of some systematic effects (energy level shifts) 
is inverted when the direction of the Raman beams are exchanged: for each 
up or down trajectory, the Raman beams directions are reversed and we 
record two velocity spectra. Finally, each determination of h/m^ and a is 
obtained from 4 velocity spectra. Fig. 3 shows two velocity spectra for the 
up and down trajectories. 

The determinations of h/m^ and a have been derived from 72 ex- 
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perimental data point taken during four days. In these measurements, the 
number of Bloch oscillations were N u p = 430 and N down = 460. Then, the 
effective recoil number is 2(N u p + N down )=1780. The dispersion of these 
n~72 measurements is x 2 /{n — 1) = 1.3 and the resulting statistical rel- 
ative uncertainty on h/m^h is 8.8 x 10~ 9 . This corresponds to a relative 
statistical uncertainty on a of 4.4 x 10 -9 . All systematic effects affecting 
the experimental measurement have been analyzed in detail in reference. 19 
The total correction due to the systematic effects is (10.98 ± 10.0) x 10~ 9 
on the determination of /i/mRb- With this correction, we obtain for a: 

or 1 = 137.035 998 84 (91) [6.7 x 10~ 9 ] (7) 

This value of the fine structure constant is labeled /i/m(Rb)2005 on Fig. 1. 




Fig. 4. Scheme of the interferometer used for the measurement of h/m^. The first pair 
of 7r/2 pulses produces a fringe pattern in the velocity distribution which is measured by 
the second pair of 7r/2 pulses. Between these two pairs of pulses, the atoms are accelerated 
upwards or downwards. The solid line corresponds to the atom in the F = 2 state, and 
the dashed line to the F = 1 state. 



4. Measurement of the fine structure constant by atom 
interferometry 

We describe in this section the results obtained in the 7r/2— 7t/2— BO— ir/2 — 
7r/2 configuration. The scheme of this intcrfcromctric method is shown on 
Fig. 4. The frequency resolution is now determined by the time Tr within 
each pair of pulses while the duration of each 7r/2 pulse determines the 
spectral width of the pulses and the number of atoms which contribute 
to the signal. This interferometer is similar to the one of the reference, 9 
except that effective Ramsey fc-wavevectors point in the same direction. 
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Consequently, this interferometer is not sensitive to the recoil energy, but 
only to the velocity variation due to Bloch oscillations which take place 
between the two sets of tt/2 pulses. 
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Fig. 5. Velocity spectra obtained in the tt/2 — tt/2 — BO — tt/2 — tt/2 configuration. 
Here TVi and N2 are respectively the number of atoms in F = 1 and F = 2 after the 
acceleration process. The spectrum on the left corresponds to the downwards acceleration 
(600 Bloch oscillations) and on the right to the up acceleration (400 Bloch oscillations). 
The frequency difference between these spectra corresponds 2000 recoil velocities. 



As in the tt — BO — tt configuration, a value of h/mnb is deduced from 
four spectra obtained with the upwards or downwards acceleration and by 
exchanging the directions of the Raman beams. An example of two spectra 
is shown on Fig. 5. In this case, the total number of Bloch oscillations is 
N u P + N down _ 10 qq ; corresponding to 2000 recoil velocities between the 
up and down trajectories. The duration of each tt/2 pulse is 400/is and the 
time Tr is 2.6 ms (the total time of a pair of tt/2 pulses is 3.4 ms). For 
these experiments, the blue detuning of the Raman lasers is 310 GHz. By 
comparison with the tt — BO — tt configuration (see Fig. 3) the resolution is 
better: the half period of the fringes is about 160 Hz when the line width of 
the spectra of Fig. 3 was about 500 Hz. Nevertheless, the reduction of the 
uncertainties is not in the same ratio. This is due to the phase noise of the 
Raman laser which becomes more important. To lower this effect, we have 
set up an active anti vibration system: then the precision of each frequency 
measurement increases by about a factor of two. After this amelioration, 
to improve the resolution, it is tempting to increase the time Tr between 
the tt/2 pulses. Nevertheless, another limitation appears, which is the size 
of the vacuum cell. Indeed, when we make the selection, the velocity of 
the atoms is close to 2Nv r and, during the selection, the atom travels a 
distance of 2Nv r Tft. For example, in the case of the spectra of Fig. 5, the 
total displacement of the atoms is about 88 mm and 53 mm for the down 
and up trajectories (the upper window of the vacuum cell is 70 mm above 
the center of the cell): practically, with the parameters corresponding to 
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the spectra of Fig. 5, we use all the size of the cell. Consequently, if we 
want to increase the delay Tr, we have to reduce the number N of Bloch 
oscillations and there is no benefit. 

To surpass this limit, we have developed a method, called atomic ele- 
vator, to better use the volume of the vacuum cell. The idea is to move 
the atoms to the top or the bottom of the cell before making the sequence 
described above. Then, we use the total size of the cell to accelerate and 
decelerate the atoms. To displace the atoms, we accelerate the atoms with 
300 Bloch oscillations during 4 ms and, after a dead time of 13 ms, we 
decelerate the atoms with 300 Bloch oscillations. This sequence displaces 
the atoms by about 50 mm. With this technique, we have increased at the 
same time N to 800 and Tr to 5.7 ms. Fig. 6 shows two records obtained 
with these parameters. The visibility of the fringes is similar to the one of 
Fig. 5 and the half period of the fringe is about 90 Hz. Now the frequency 
difference between the two spectra corresponds to 3200i; r . During the selec- 
tion, the atomic velocity is about 10 m/s and the atom travels a distance 
close to 6 cm. This shows the limitation due to the size of the vacuum cell. 

6 = 23 819 247.9 (0.5) Hz 5 = -24 473 138.0(0.6) Hz 




' 0.15 

100 200 300 400 600 700 800 900 1000 

5 meas - 8 s e/ - 23 819 000 (Hz) S meas -5 se , + 24 474 000 (Hz) 



Fig. 6. Velocity spectra obtained in the iz/2— iz/2— BO— 7r/2— tt/2 configuration with the 
atomic elevator. The spectrum on the left corresponds to the downwards acceleration (800 
Bloch oscillations) and on the right to the upwards acceleration (800 Bloch oscillations). 
The frequency difference between these spectra corresponds 3200 recoil velocities. 



We present now the result deduced from 221 measurements of h/m^va- 
For these measurements, we have used the two methods described previ- 
ously, with and without the atomic elevator. The total number of Bloch 
oscillations N u p + N down varies from 200 to 1600. The dispersion of these 
n = 221 measurements is % 2 /(n — 1) = 1.85 and the resulting relative sta- 
tistical uncertainty on a is 3 x 10~ 9 . The systematic effects arc similar to 
the ones described in reference. 19 The two main effects are due to the geom- 
etry of the laser beams and to the second order Zeeman effect. To evaluate 
the first effect, we have measured the wave front curvatures with a Shack- 
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Hartmann wave front analyzer (HASO 128 from Imagine Optics). From 
these measurements, we have obtained a correction of (—11.9 ± 2.5) x 1CP 9 
on the determination of a. As explained above, the effect of parasitic level 
shifts is eliminated by inverting the direction of the Raman beams. Never- 
theless, this assumes that the measurements are made exactly at the same 
position when the direction of the Raman beam is inverted. In fact, these 
positions are not exactly the same because the directions of the recoils given 
at the first Raman transition are opposite. For the timing used in our ex- 
periment, they differ by about S z = 300 ^m. We have precisely measured 
the spatial magnetic field variations to control this effect. This correction 
depends of the number of Bloch oscillations. For example, in the case of the 
records of Fig. 6, its value is (7± 1) x 10~ 9 . Finally the relative uncertainty 
due to the systematic effects is 3.4 x 10~ 9 and we obtain for a: 

a' 1 = 137.035 999 45(62) [4.5 x 10~ 9 ] (8) 

This value corresponds to the point labeled ft,/m(Rb)2008 on Fig. 1 and 
is in agreement with our 2005 measurement. Our two results are also in 
agreement with the most precise value deduced from the electron anomaly 
(labeled a e (Harvard) on Fig. 1). 

5. Conclusion 

We have presented two determinations of the fine structure constant a. De- 
pending on the Raman pulses arrangement (tt — BO — tt or tt/2 — n/2 — 
BO — 7r/2 — 7r/2 configurations), our experiment can run as an atom in- 
terferometer or not. The comparison of the two resulting values, which are 
in good agreement, provides an accurate test of these methods. The com- 
parison with the value extracted from the electron anomaly experiment 6 
is either a strong test of QED calculations or, assuming these calculations 
exact, it gives a limit to test a possible internal structure of the electron. 
Our goal is now to reduce the relative uncertainty of a. We are building 
a new experimental setup with a larger vacuum chamber. With the new 
cell, we plan to multiply the number of Bloch oscillations by a factor of 
three. Then, it will be possible to reduce the uncertainty at the 10~ 9 level 
to obtain an unprcccdcnt test of the QED calculations. 
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